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ABSTRACT 

We present Stellar Locus Regression (SLR), a method of directly adjusting the instrumental broad- 
band optical colors of stars to bring them into accord with a universal stellar color-color locus, produc- 
ing accurately calibrated colors for both stars and galaxies. This is achieved without first establishing 
individual zeropoints for each passband, and can be performed in real-time at the telescope. We 
demonstrate how SLR naturally makes one wholesale correction for differences in instrumental re- 
sponse, for atmospheric transparency, for atmospheric extinction, and for Galactic extinction. We 
perform an example SLR treatment of SDSS data over a wide range of Galactic dust values and 
independently recover the direction and magnitude of the canonical Galactic reddening vector with 
14-18 mmag RMS uncertainties. We then isolate the effect of atmospheric extinction, showing that 
SLR accounts for this and returns precise colors over a wide of airmass, with 5-14 mmag RMS residu- 
als. We demonstrate that SLR-corrected colors are sufficiently accurate to allow photometric redshift 
estimates for galaxy clusters (using red sequence galaxies) with an uncertainty a{z)/{l + z) = 0.6% 
per cluster for redshifts 0.09 < z < 0.25. Finally, we identify our objects in the 2MASS all-sky 
catalog, and produce «-band zeropoints typically accurate to 18 mmag using only SLR. We offer open- 
source access to our IDL routines , validated and verified for the implementation of this technique, at 
|http://stellar- locus-regression. googlecode.com. 

Subject headings: galaxies: fundamental parameters — methods: data analysis — stars: fundamental 
parameters — stars: statistics — techniques: photometric 



1. INTRODUCTION 

The observed broadband colors — i.e. fiux ratios — of ce- 
lestial objects depend on the photon spectral energy dis- 
tribution of the source; on extragalactic. Galactic and 
atmospheric scattering and absorption along the line of 
sight; and on the instrumental sensitivity function over 
the wavelengths of interest. One challenge of astro- 
nomical photometric analysis is to disentangle, from a 
given set of observations, the source's colors from such a 
plethora of perturbing factors. 

Colors hold information about a source's type, tem- 
perature, metallicity, and redshift. A source's apparent 
magnitude, on the other hand, also depends on its dis- 
tance and on the size and nature of the emitting regions. 
We assert that for almost all astrophysical endeavors, ac- 
curate photometric colors are more useful than high ac- 
curacy magnitudes, especially because we seldom know 
distances well enough to convert from apparent to abso- 
lute magnitudes at the percent level. Furthermore, for 
photometric redshift techniques that use a prior on mag- 
nitudes, the broad luminosity function of galaxies gener- 
ates a span in magnitude that far exceeds the range in 
color, for a given galaxy type at a particular redshift. 

The current and next generation of wide-field 
multicolor survey projects include the CFHTLS"'^, 
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PanSTARRS^, BCS^, DES^ LSST^, and SkyMapper^ 
Object classification and distance estimation with pho- 
tometric redshifts are necessary starting points for ex- 
tragalactic science using these surveys, and these in turn 
depend on knowing calibrated colors. The technique we 
describe here can be used very early in a survey to obtain 
highly accurate colors of objects, as well as magnitude 
estimates. This should allow more rapid exploitation of 
new survey data — indeed, of nearly any multiband data. 

The standard approach to determining colors of 
sources is to first calibrate magnitudes in all observed 
passbands (such as 5, r, i, and z), and then subtract the 
calibrated magnitudes to obtain calibrated colors {g — r, 
g — i, etc). This is typically time-consuming, both at the 
telescope and in the analysis phase. Establishing pho- 
tometric zeropoints for a stack of multiband images re- 
quires separate observations of spectrophotometric stan- 
dard stars to measure the instrumental sensitivity and 
estimate atmospheric extinction, and spatiotcmporal in- 
terpolation of the calibration parameters to the science 
fields under "photometric conditions." 

In this paper we describe how to calibrate colors di- 
rectly from objects cataloged from multiband, flat-fielded 
images of a field, without having to first determine the 
corresponding photometric zeropoints, and without the 

^ http: //pan- Starrs . if a.hawaii . e du/public/] 
^ http://cosmology.uiuc.edu/BCS/ 
^ https : //www. darkenergysurvey . org/ 
^ http : // www . Isst . org/Science/lss t_baseline ■ shtml| 
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usual repeated measurements of standard stars. We 
demonstrate how this can be done accurately, yielding 
colors accurate to a few percent, and rapidly, allowing 
for optimal use of allocated telescope time. 

Our technique exploit s the optical and infrared color- 
color stellar locus (cf. ICovev et al.l l2007l : llvezic et"al] 
l200l . the one-dimensional and astrophysically funda- 
mental track that stars occupy in color-color space. The 
majority of stars lie somewhere along this locus, at a 
position that depends primarily upon effective tempera- 
ture. The universality of the stellar locus was exploited 
in the Oxford- Dartmouth Thirty-Degr ee Survey ( "stellar 
locus fitting," iMacDonald et all l2004l ) to stabilize pho- 
tometric zeropoints in non-photometric conditions, and 
in Sloan Digital Sky Survey 's (SDSS) Stripe 82 ("stellar 
locus method," llvezic et al.1 ,2007') to account for differ- 
ences in the response function of different detectors in 
the SDSS instrument. Our approach is different from 
these mainly because we do not first establish photomet- 
ric zeropoints per band: we immediately calibrate colors 
in all fields, and only optionally solve for the calibrated 
apparent magnitudes, using the stellar locus. We do not 
use SLR as a diagnostic or corrective tool, but as our 
primary calibrator. 

Our Stellar Locus Regression (SLR) approach builds 
on these previous pioneering works. We establish the 
location of the stellar feature in instrumental color- 
color space, and we determine what transformations are 
needed to bring this into coincidence with the known 
location of a standard stellar locus. Applying appro- 
priate color-corrections to the entire rest of the catalog 
automatically accounts for all of the standard calibra- 
tion terms, including zeropoints, atmospheric extinction, 
aperture corrections, and Galactic extinction. The tech- 
nique is straightforward, fast, and allows observers to 
forgo the usual standard star observations altogether be- 
cause all observed stars are expected to lie along the same 
stellar locus. We have calibrated fields with SLR using 
as few as 7 stars in fields of view as small as 4' x 8'. 

Figure [T] schematically illustrates the technique. We 
perform SLR on new data (see ^4.2114^ taken with 
the IMACS instrument ([Dressier et al.li2003t lOsip et all 
[2OO81) on the Magellan 6.5 m telescope. All panels show 
our adopted standard stellar l ocus line and stella r den- 
sity contours, reproduced from' Covev et al] (|2007l . W2.1\i . 
along with 36 IMACS stellar colors (red points). The top 
panels show instrumental stellar colors with the stan- 
dard locus. We perform SLR, neglecting instrumental 
color terms, with results shown in the middle panels. 
We perform SLR again after measuring color terms in- 
dependently from an external standard star field, with 
results shown in the bottom panels. By way of illus- 
tration, the vectors in the middle panels show the ex- 
pected direction and magnitude of extinction by Galac- 
tic dust (Ar = 0.2 mag and Ry = 3.1, estimated using 
ISchlegel et al.| [l998) and the atmosphere (1.3 airmasses). 

We adopt a standard stellar locus in fj^l and then mo- 
tivate and build a picture of what comprises it so as to 
understand its universality. In ^J3]we outline the math- 
ematics of color calibration and describe an optimal al- 
gorithm for real-time color calibration, which we apply 
to real data in a series of tests O. In H4.ll we per- 
form our first test, applying SLR to already-calibrated 
SDSS photometry, where we recover the canonical Galac- 
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Fig. 1. — An illustration of Stellar Locus Regression (SLR). Col- 
ors are plotted on the SDSS photometric system. All panels show 
the standard stell ar locus ( black line and gray density contours), 
reproduced from C ove y et al . (SOOT!) . Red points are stellar colors 
obtained from a Source Extractor analysis of flat-fielded Magellan 
6.5 m IMACS images. Top panels: The instrumental IMACS colors 
are plotted, with a clear mismatch between them and the standard 
locus. Middle panels: SLR is performed with only a common trans- 
lation vector applied to the instrumental colors. Note the color- 
dependent discrepancies in the upper right portions of the central 
panels. By way of example, the vectors show the expected direc- 
tion and magnitude of extinction by Galactic dust {Ar = 0.2) and 
the atmosphere (1.3 airmasses). Bottom panels: Color terms are 
measured from a single observation of a field containing standard 
stars. Fixing these color terms, a new best-fit translation is deter- 
mined, which brings the observed colors onto the SDSS-calibrated 
color system, as defined by the stellar locus. This SLR analysis, 
when the corrections are then applied to all objects in the photo- 
metric catalog, allows us to rapidly obtain highly accurate colors 
on the SDSS system, directly from fiat-fielded data, with a single 
correction step that accounts for atmospheric extinction, Galactic 
extinction and instrumental response differences. 

tic reddening vector in direction and magnitude over a 
wide rage of dust thicknesses. Section 14.21 examines the 
fundamental repeatability of SLR color and magnitude 
measurements by comparing SDSS data to data from a 
different instrument. Section 14.31 isolates the effect of 
atmospheric extinction. Our final test appears in ^14. 41 
where we recover the spectroscopic redshifts of 11 low 
redshift galaxy clusters to high accuracy using only SLR 
colors. We end with a discussion (SJ5]) and conclusions 

m- 

2. THE STELLAR LOCUS 

The nearly-blackbody emission spectra of stars place 
them predominantly along a line in optical and infrared 
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color-color space. Stellar color therefore depends pri- 
marily on effective temperature and is the basis for 
the Morgan-Keenan (MK) stellar classification system 
(|Morgan et al.|[r943f) . Real stellar atmospheres deviate 
from blackbody behavior because of molecular absorp- 
tion and metallicity effects. Observations, which capture 
stellar light that has traversed the dust of our galaxy, 
our atmosphere, and our telescope systems, will produce 
stellar color loci that deviate further from the expected 
intrinsic behavior. We have based our SLR approach on 
the assumption that we can identify a stellar locus in 
grizJHK color space that is intrinsically universal. We 
explore the extent to which this is true by inspecting 
what comprises an observed stellar locus. 

2.1. A Standard Locus 

We adopt as our st andard the empirical stellar locus 
of ICovev et al.l (|2007f l. Those authors calculated the 
running-median of the colors of ^ 10^ stars, from high 
quality observations by both the SDSS and 2MASS sur- 
veys. The line-of-sight Gal actic dust for their sa mple was 
estimated from the maps of lSchlegel et"an (|1998l . SFD) to 
be Ar < 0.2. We additionally smooth this locus line with 
a 0.2 mag size boxcar averaging kernel to reduce some of 
the scatter between adjacent points. These are the data 
we present in Figure [T] — we call this the standard stellar 
locus line. 

The two most salient features of the standard locus 
are its nearly one-dimensional nature and a dramatic 
break or kink in the (g — r,r — i) plane a,t r — i 0.7. 
As we explore below (i i2.2[) . blue- ward of the kink are 
mainly evolved and main sequence (MS) A- through 
K-type stars, and re d -ward are ma i nly M-type dwarfs 
(IFinlator etal] 120001: jHawlev et all 120021 : ICovev et all 
l2nn7H,Turic et al.ll200af) . 

Because of the well-defined relationships between the 
color (effective temperature), age, and intrinsic luminos- 
ity of stars, each point along the stellar locus probes 
roughly predictable 3-dimensional spatial volumes, given 
some fixed dynamic range for the observations. The red- 
dest MS stars that make up the r — i> 0.7 branch of the 
stellar locus are less luminous, so the effective volume an 
observation probes is small and nearby. Likewise, bluer 
MS and evolved stars are intrinsically brighter, so the 
effective volumes probed at those colors are larger and 
farther away — and because both MS and evolved stars 
constitute the r — i < 0.7 portion of the stellar locus, 
observations are sensitive to a plurality of volumes and 
distances. 

A comprehensive Galactic structure and population 
synthesis analysis is beyond the scope of this paper, but 
we will explore below the factors that give rise to unifor- 
mity in the observed stellar locus, including the dramatic 
color-volume effects at typical SDSS depths. 

2.2. Color- Volume Effects 

Each pointing of an astronomical camera images a 
cone-shaped region of space, set by the solid-angle field 
of view O of the instrument. For a single frame there is 
also a limited observable dynamic range in apparent mag- 
nitude: the brightest observable objects are determined 
by the saturation limit of the system, and the faintest 
useful objects must satisfy some selection in signal-to- 
noise ratio. A star of a particular absolute magnitude 



TABLE 1 

Detectability Volumes vs. Absolute Magnitude. 



Mr'' -Rinner*" -Router'^ Detectability Volumc'' 
(AB mag) (pc) (pc) (arb units) 

-5 
-4 
-3 
-2 
-1 

1 
2 
3 
4 
5 
6 
7 
8 
9 
10 
11 
12 
13 
14 
15 



^ Absolute magnitudes of stars in the SDSS r- 
band.'' Closest distances (in parsecs) for the object 
to be observable, for a notional saturation limit of 
r = 14. Furthest distances (in parsecs) for a lOcr de- 
tection limit of r = 22. Volume (in arbitrary units) 
within which objects of a given absolute magnitude 
can be detected, subject to the instrumental dynamic 
range constraints. 

and corresponding color is therefore detectable within a 
truncated cone of opening angle f2, with an near-edge de- 
termined by the saturation limit and a far-edge set by the 
required signal-to-noise ratio. For simplicity we assume 
an unextincted line of sight. 

Table [T] illustrates the dramatic selection-effect of dy- 
namic range on observed stars of various intrinsic lumi- 
nosities. We take a notional saturation limit of r — 14 
and a faint detection limit at r = 22, which are approxi- 
mate values appropriate for SDSS. For a dynamic range 
of 8 mag, the outer detection edge of the detectability 
cone is always ^ 39 times farther away than its inner 
saturation edge, for a given absolute magnitude object. 
Table [T] also shows that for each additional magnitude in- 
crease of stellar luminosity, the survey volume increases 
by a factor of four. Luminous stars are detectable over 
a vastly larger volume than fainter stars, since the outer 
edge of the detectability region is proportional to the dis- 
tance to the inner edge. 

As a concrete example, at SDSS depths and Galactic 
latitudes |6| = 90 deg, absolute magnitudes M^ > 12 are 
detectable only closer than ^ 1 kpc, wh ich is the height 
of the Galactic disk (jJuric et al.l I2008D . Objects with 
M,. < 4 are seen only at Galactic heights \Z\> Ikpc, so 
are in the halo. 

Absolute stellar magnitude varies with stellar color 
and age in well understood ways. To explore the 
relation in the SDSS magnitud e system, w e used 
the model ste llar populations of 'Marigo et al.| (|2008l 
Ihttp : //stev . oapd . inaf . it/cgi-bin/ cmd_2 .1). This 
online tool modeled stars with initial masses down to 
0.15 Af©. We have just shown that the typical SDSS dy- 
namic range is sensitive to both halo and disk stars, so 
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Fig. 2. — Comparisons of modeled Mr vs. g — r color-magnitude 
diagrams (left column), r — i vs. g — r (central column), and i — 
z vs. g — r (rightmost column) color-color diagrams for different 
stellar populations, for zero Galactic extinction. The upper row 
corresponds to a stellar halo population with [Fe/H] = —1.5 and 
an age of 13 Gyrs. The central row also corresponds to a 13 Gyr 
age, but with a characteristic disk metallicity; [Fe/H] = —0.5. The 
bottom row is what one might expect for more recent star formation 
in the Galactic disk, with [Fe/H] = —0.5 and an age of 5 Gyr. 
Absolute magnitudes are indicated by the angle of the whisker, 
with 9 = 9Mr mag~^ deg. Objects with r — i > 0.7 are shown 
as triangles. Only stellar populations with [Fe/H] > —0.7 produce 
the excursion to r — i > 0.7, seen in the central and lower rows, 
that is essential to the SLR method. 



we consider three models, using llvezic et"aLl (|2008() as a 
guide: 

1. a 13 Gyr lialo population, with a metallicity 

[Fe/H]haio = -1.5, 

2. a 13 Gyr disk population of higher metallicity, 

[Fe/H]disk = -0.5, and 

3. a 5 Gyr disk population of higher metallicity, 
[Fe/H]disk - -0.5. 

Figure [2] shows the color-magnitude and color-color dia- 
grams that illustrate the connection between absolute 
magnitudes and regions in color-color space for these 
populations. 

Figure [2] shows that only disk-metallicity stars popu- 
late the locus red-ward of the kink. The main sequence 
M dwarfs in this color range have Mr ^11 mag and so 
in the SDSS dynamic range these objects lie at distances 
< 1.5 kpc from Earth. These stars are almost entirely 
in the disk. The IVI giants red-ward of the kink have 
Mr < —2 mag corresponding to distances > 15 k pc. Few 
disk-metallicity stars lie at such distances (Ivezi c et al.l 
|2008| ). Therefore the vast majority of the stars that pop- 
ulate the (r — i) > 0.7 region of t he loc us at SDSS depths 
are disk dwarfs. iHawlev et al.l ()2002f ) also support this 



viewpoint, attributing SDSS objects with 1 < (r — i) < 2 
to JVIO through JVI6 stars at typical distances ranging from 
1 to 0.1 kpc respectively. 

Continuing these arguments, stellar locus stars blue- 
ward of the {r — i) ^ 0.7 kink comprise of a blend of 
JVIS disk stars, JVIS halo stars, and evolved (post-turnoff) 
halo stars. The luminosity of the bluer MS stars are 
4 mag < Mr < 11 mag and are detectable from both 
the disk and halo. We also see that even for a stellar 
population 13 Gyrs old, with a main sequence turnoff at 
around Air ~ 4, blue turnoff stars that are closer than 
1 kpc would saturate the detector. 

The relative proportion of disk and halo stars that 
make up the stellar locus at (r — i) < 0.7 (blue- ward 
of the kink) is magnitude dependent, llvezic et al.l ()2008D 
show that for 0.2 < {g - r) < 0.4, at g ~ 17.5 half the 
stars have metallicity characteristic of the disk, and half 
are halo stars. At fainter magnitudes the proportion of 
halo stars increases. 

Armed with this understanding of the observed stel- 
lar locus's constituent populations, we explore the effect 
their corresponding metallicity distributions should have. 

2.3. Metallicity 

A number of previous authors have explored the influ- 
ence of metallicit y on stellar color s in the SPSS pho- 
tometric sy stem (iLenz et al l Il998t Covev et alJ 120071 : 
lJuric et al.ll2008l : llvezic et al.fl2008D . llvezic et al.r(|2008f ) 



demonstrated that most halo stars have a common metal- 
licity that is reasonably well described by a Gaussian 
distribution with median [Fe/H]haio = —1.5 and width 
Chaio = 0.3 dex. The median disk metallicity was shown 
to be about —0.5 and exhibits half the scatter of the halo, 
with CTdisk = 0.16 dex. 

The biggest observable broadband effect from metal- 
licity occurs in the u-band, which is one of the reasons 
we exclude that band from consideration here. Stellar 
metallicity effects in g — r are ten ti mes smaller than in 
u—g. Using this fact and Figure 21 of llvezic et al.l |2008), 
the expected spread in g — i color for halo stars should 
be of order 0.05 mag. This is a conservative upper limit 
estimate of the intrinsic width of the halo contribution 
to the stellar locus from metallicity alone. 

We have modeled the effect of metallicity on the stel- 
lar MS, and present some of the results in Table |2l 
We used a grid of model atmospheres deriv ed from 
the Phoenix project (jBrott fc Hauschildtll2005l ) and the 
iKurucd ()19 93) library for the typical disk and halo metal- 
hcities, [Fe/H] = —0.5 and —1.5 dex. The range of pa- 
rameters used was 3, 500X < T^s < lO.OOOiiT for M5 
to AO using the Phoenix models and 15. OOOK < Tcs < 
45,000A' for B5 to 05 from the iKuruczl (jl993h models, 
and we adopted a typical main sequence log(g) for each 
spectral type. We then calculated the broadband colors 
based on the SDSS filter and CCD throughput. 

We found that the greatest influence of metallicity on 
broadband color is in the late-type stars (around MO and 
later) and is in opposite direction to the color displace- 
ment of early-type stars. The relatively cool atmospheres 
of the red stars are rich with molecules, and the molec- 
ular absorption features blend together with other lines 
to form a continuum that differs substantially from a 
blackbody spectrum. We also observed that a minimum 
metallicity of [Fe/H] = —0.7 is needed for red stars to 
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TABLE 2 
Model Stellar Color Difference 
Between Z=— 0.5 and Z=— 1.5. 



MK type 


5(g - 


(S(r — lY 


6{i - 2)=- 


05V 


-0.005 


-0.002 


-0.003 


BOV 


-0.000 


-0.002 


-0.002 


B5V 


-0.003 


-0.002 


-0.005 


AOV 


-0.000 


-0.004 


-0.003 


A5V 


0.004 


-0.007 


-0.007 


FOV 


0.004 


-0.010 


-0.012 


F5V 


0.019 


-0.009 


-0.013 


GOV 


0.031 


-0.012 


-0.015 


G5V 


0.037 


-0.013 


-0.015 


KOV 


0.053 


-0.015 


-0.014 


K5V 


0.124 


0.002 


-0.014 


MOV 


0.004 


0.067 


0.031 


M5V 


-0.137 


0.117 


0.080 



^ As a comparison, for galactic extinction 
{Ar = 0.2) the color perturbation is (0.076, 
0.048, 0.044); for 1.3 atmospheres of extinc- 
tion, the color perturbation is (0.103, 0.053, 
0.010). 



populate the locus red-ward of the kink. Less metallic- 
ity "irons out" the kink, making the model colors col- 
inear with the bluer stars. This observation is consis- 
tent with the conclusions drawn in W2.2\ A continuum 
of metallicities near this threshold would smear the ob- 
served stellar locus in the redder region. The apparent 
lack of such smearing in the observed locus indicates the 
pa ucity of stars at this intermediate metallicity, as found 
by llvezic et all (|2008[ ). It is therefore the high mean 
metallicity and molecular absorption of disk stars that 
produce the kink in the SDSS stellar locus in the first 
place, and that make our implementation of SLR possi- 
ble. 

2.4. Age 

Younger stellar populations have a brighter, bluer main 
sequence turnoff point, and this shifts the blue tip of the 
stellar locus towards the blue. This does not distort the 
basic shape of the stellar locus. For very young popula- 
tions at high metallicity there is also a red-ward exten- 
sion of the main locus, due to massive stars evolving to 
r — i > 0.7, but we concluded in ^ 32. 21 that these stars 
do not appear in the stellar locus at SDSS depths, since 
these young stars would be in the disk and would sat- 
urate. Other features in the color-magnitude diagram, 
such as helium-burning loops, essentially remain confined 
to the standard locus line and do not change its shape. 
Our methods (Sj3|) are therefore expected to be largely 
insensitive to population age effects. 

2.5. Binaries and Unresolved Superpositions 

Superpositions of stars, either as physical binaries or 
chance coincidences, most often fall within the standard 
locus area, but can also produce systematic outliers and 
therefore change the observed stellar locus shape. For 
the latter to happen, the objects within a single PSF 
must have different colors, but similar magnitudes. One 
can adopt two philosophical resolutions to this issue. 

One approach is to realize that both our standard locus 
and the instrumental colors should contain binaries, and 
so the effect will average to zero over many fields. While 
this may be true statistically, if the stars drawn from 



an observed field are few in number, the observed stellar 
locus may be distorted by small number statistics. 

Another approach is to iteratively excise stellar colors 
that are significantly different from the standard color 
expectation. This should suppress locus shape distor- 
tions due to the occasional unresolved binary, but nev- 
ertheless thejncidence of bin ary outliers is expected to 
be low. ISmolcic et al.l (|2004f ) have shown the existence 
of an "echo" of the stellar locus, which they ascribed 
to unresolved binaries of similar luminosity but different 
effective temperatures. They assessed the number of ob- 
jects in this "echo locus" as being fewer than 1/2000 as 
numerous as the objects that occupy the main stellar lo- 
cus. The combined low incidence and known location of 
these outliers makes these binaries unproblematic. 

2.6. Variable Stars 

Photometric variability is a source of systematic er- 
ror for stellar locus color methods. If a star's observed 
brightness varies in the time between observations in the 
various passbands under analysis, even if the underly- 
ing stellar color does not vary, the magnitude differences 
will introduce a shift in color-color space. Since the SDSS 
photometric data in different bands are obtained within a 
few minutes of each other, the SDSS-band standard stel- 
lar locus should have minimal contamination from this 
effect. Similarly, the 2MASS instrument obtained J, H, 
and K photometry essentially simultaneously. Combin- 
ing 2MASS and SDSS photometry, however, can fall prey 
to variability, since the observations were taken at differ- 
ent times. So establishing the joint optical/NIR standard 
stellar locus must attend to this issue. 

Similarly, if the images under analysis are obtained at 
different times, variability on this ti mescale will distort 
the derived calibration color shifts. ISesar et al.l ()2007l ) 
explored the incidence of variables in the SDSS color- 
color diagram. They showed that there is a substantial 
variation in the fraction of objects that exhibit variabil- 
ity, across the color-color diagram. The blue tip of the 
stellar locus has two classes of variables: low redshift 
QSOs and RR Lyrae stars. These can be suppressed 
by selecting a judicious region in color-color space where 
the standard locus is matched to the instrumental color 
distribution, iteratively. 

The mul t i-epoc h analysis of SDSS photometry in 
ISesar et al.l (120071 ) indicates that fewer than 5% of the 
objects in the main stellar locus exhibit variability with 
Ag > 0.05 mag. A sigma-clipped iterative analysis of 
multi-epoch photometry should be able to produce a 
cleansed standard stellar locus. 

2.7. Galactic Dust 

Galactic dust is another extrasolar sou rce of stellar lo- 
cus pe rtur bations. The dust is shown by [Marshall et al.l 
(120061) and lJuric et all (|2008f l to be confined to a sheet 
roughly 100 pc above and below the Galactic midplane. 
At SDSS depths, the vast majority of stars are therefore 
behind the dust. 

To first approximation the extinction obeys the canon- 
ical Rv = 3.1 reddening law and the degree of extinction 
Av follows the maps of SFD. This induces a simple over- 
all color-color vector shift — reddening — whose direction 
depends on Ry a nd whose magnitude depends on Ay- 
llvezic et al.l (|2007[ ) showed that the stellar locus position 
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reflects these predictions in some regions of SDSS Stripe 
82, but breaks down as Galactic latitude |5| decreases 

(S3)- 

The usual adoption of extinction coefficients from the 
SFD appendix assumes an underlying spectral energy 
distribution typical of an e lliptical galax y. This is not 
valid for Galactic stars fee lMcCallll2004D . To next ap- 
proximation, a more correct treatment would provide ex- 
tinction coefffcicnts for each stellar spectral type. Errors 
of this form induce differential distortions of the stellar 
locus as a function of effective temperature. 

Other deviations occur when the basic assumptions of 
extinction behavior break down. If the redde ning law is 
not strictly Ry = 3.1 fee iLarson" fc Whittetl 2005), the 
stellar locus shifts in different directions. If some ob- 
served stars are in front of the dust and others behind, 
the stellar locus will show additional scatter, preferen- 
tially in the reddening directions. Because the stars are 
most often behind the dust. Galactic extinction is a ma- 
jor contributor to the observed stellar locus properties. 

The SLR approach makes no specific assumption about 
any value of Ry- It corrects for Galactic extinction sub- 
ject to the assumption that the observed stars all lie be- 
hind a common Galactic extinction layer. We suspect a 
fruitful approach might be to make an initial adjustment 
to the instrumental colors using SFD and Ry = 3.1, and 
then run SLR on the resulting catalog. 

2.8. Summary 

Taking all this information together, we have clear ex- 
pectations for what constitutes the typical observed stel- 
lar locus, and what factors change its observed position 
or shape. 

Stars with r — i > 0.7 arc mainly faint M dwarfs in the 
disk, with correspondingly high metallicity. Since these 
objects are only visible out to ~ 1 kpc, only metallicity 
variations in this local region of the Galaxy could perturb 
the kink region of the observed locus. 

The stars blue- ward of the kink, however, are a 
magnitude-dependent combination of halo stars and disk 
main sequence stars. The metallicity dependence of this 
region of the locus is small. 

While the location of the bluest terminating edge of the 
stellar locus is an indicator of the age of the underlying 
population, it does not distort the basic shape of the 
locus line. 

Galactic dust should be a prime source for locus shifts 
and possibly small distortions. At high Galactic lati- 
tudes the vast majority of SDSS stars reside behind the 
dust lanes of the Milky Way. At low Galactic latitudes, 
even distant stars come closer to the Galactic plane. We 
would therefore expect to eventually see a smearing of 
the stellar locus due to stars suffering different amounts 
of extinction along the line of sight. 

Finally, atmospheric and instrumental effects will nat- 
urally perturb the stellar locus, and we model these in 
our mathematical formulation of Stellar Locus Regres- 
sion (^. Given what we have learned about the nature 
of the observed stellar locus, our conception of SLR be- 
gins with a simple set of assumptions. 

1. The standard stellar locus is representative of the 
typical stellar populations that we will observe in 
practice, and is sufficiently uniform so as to con- 



stitute a calibration standard. The standard locus 
stars: 

(a) are disk dwarfs in the red and both disk and 
halo stars in the blue, thanks to the dramatic 
color-volume effects at SDSS depths, and 

(b) lie at high enough Galactic latitude to put 
them behind the dust, but suffer low dust ex- 
tinction. 

2. Stars we observe in practice are always behind the 
Galactic dust sheet. This way, SLR directly out- 
puts dereddened colors. 

3. The Galactic dust extinction can obey any power 
law Ry. 

4. The Galactic dust extinction is locally smooth. 

5. Observations include the g-band, and are deep 
enough to observe at least a few of the faint M 
dwarfs that constitute the k ink feature (limiting 
g > ISmag. ll^ic et al.l[200i ). 

6. The input images from which instrumental pho- 
tometry is extracted are properly flat-fielded. 

Many of these are testable in isolation (21) > but some pos- 
sibly degenerate effects, such as simultaneously anoma- 
lous metallicities and extinction laws, may be difficult 
to disentangle. We deduce nonetheless that the combi- 
nation of magnitude dynamic range selection effects and 
the relative insensitivity of the grizJHK stellar locus to 
stellar metallicity, age differences, and binary and vari- 
able contamination produces an observed stellar locus 
that is uniform enough to achieve calibration of colors 
using a standard stellar locus. 

3. THE METHOD 

The basis of Stellar Locus Regression is to transform 
instrumental stellar colors so that they align with a stan- 
dard locus, on the SDSS photometric system. This re- 
quires typical data preprocessing and photometry that 
produces instrumental stellar colors from single-epoch, 
ffat-fielded images. The only calibration images required 
for SLR are a single set of multiband observations of a 
high-density standard star field, from which instrumen- 
tal color terms are measured. Periodic updates to the in- 
strumental color terms will depend on the timescale over 
which these terms evolve. Mosaic imagers will benefit 
from chip-by-chip color terms. Then, with fixed color 
terms applied, SLR calibrations are performed on sci- 
ence frames by iteratively transforming the instrumental 
stellar colors to optimize a goodness-of-fit (GOF) statis- 
tic. The resulting best-fit parameters, including the in- 
strumental color terms, define the color transformation 
that achieves the SLR calibration. Uncertainties in the 
calibration are estimated numerically. These calibration 
terms are then applied to all the cataloged photometry of 
objects appearing in the same images. This way all ob- 
jects in the field are calibrated using the same stars lying 
in that field, and no spatial nor temporal interpolation 
is required. 

In this section we describe and motivate our calibra- 
tion equations, discuss our chosen GOF and method of 
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error estimation, and outline a practical algorithm that 
we have implemented for the real-time calibration of griz 
colors. 

3.1. Color Transformations 

A color transformation is the mathematical transfor- 
mation of colors by translations, scalings, rotations and 
shears. Instrumental colors are represented with the vec- 
tor^ c, and "true" colors, which we take to be on the 
SDSS photometric system, are represented with the vec- 
tor Cq. Instrumental colors relate to the true colors via 
the color transformation equation 

c=/t+(l + B)co. (1) 

The color translation vector k accounts for first order 
atmospheric extinction. Galactic extinction, zeropoints, 
and any other additive effects, known or unknown. We 
discuss K in §3.21 The color term matrix B is populated 
by zeros and constants, typically small, and corresponds 
to the instrumental color terms from conventional pho- 
tometric calibrations. We discuss B in Appendix [Xl 
shows in mathematical detail how Equation ([T]) relates 
to photometric calibration equations that astronomers 
typically adopt. 

3.2. The Color Translation Vector 

The color translation vector k accounts for differences 
of zeropoints a, atmospheric extinction E, and Galac- 
tic extinction A. The scalar elements of the translation 
vector are 

where n and m are elements of the filter set, for example 
n,m € {g,r,i,z}, and n ^ m. See Appendix [X] for a 
derivation of this relation, starting from traditional pho- 
tometric calibration equations. 

In traditional photometric calibration, the terms on 
the right-hand side of Equation ([2]) are estimated inde- 
pendently. Atmospheric extinction is usually modeled 
using a Bouger extinction law, fc„X„, where fc„ is a filter- 
dependent constant and X„ is the airmass through which 
the i?-band image was taken. The atmospheric extinc- 
tion constants are extracted using intermittent standard 
star observations at a range of airmasses and interpo- 
lated in space and time to the science frames. This as- 
sumes that the linear airmass model reflects truth and 
is spatiotemporally invariant. Galactic extinction is nor- 
mally estimated from SFD — a procedure that assumes a 
single Rv = 3.1 dust reddening law (cf. ISchlegel et al.l 
Il998l ). Finally, zeropoints are estimated using the same 
standard star frames used to measure the atmospheric 
terms. The zeropoints standardize the photometry, and 
account for differences in instrumental throughput be- 
tween facilities and any other unmodeled additive effects. 
Zeropoints are then interpolated in space and time to the 
science frames in "photometric conditions," which again 
assumes spatiotemporal invariance. 

SLR is fundamentally different and is not subject to 
these assumptions. SLR fits for each element of the vec- 
tor K directly, and therefore calibrates the entire right- 
hand side of Equation ^ in one step. The atmospheric 

Our convention is to assign vectors boldface, lowercase letters, 
and matrices boldface, capital letters. The elements of vectors, 
matrices, and tensors, like all scalars, are not boldface. 



and dust extinctions and the zeropoints are not estimated 
independently. This is made possible by the universal- 
ity of colors of abundant MS stars (^. Because of this 
unique ability of SLR, it does not matter what mathe- 
matical form the extinction and zeropoint terms take in 
reality, as long as they are additive, because all additive 
systematic effects are accounted for together by k during 
SLR. For example, the Bouger atmospheric extinction 
law is likely incorrect in the z-band, where water absorp- 
tion is saturated and the additive airmass dependence 

1 /2 

probably better fo llows an Xz airmass power law (cf. 
IStubbs et ani2007f ). Further, as we have mentioned it 
is probable that not all Galactic dust columns obey the 
canonical Ry — 3.1 extinction law. SLR is immune to 
additive mis-modeling of the atmosphere, Galactic dust, 
overall instrument sensitivity differences, aperture cor- 
rections, and so on. 

Figure [T] demonstrates the effect of leaving the color 
translation vector free during Stellar Locus Regression. 
The top panels show how the known, standard locus 
compares to an instrumental stellar locus. In the mid- 
dle panel, we perform SLR, letting k be free and fixing 
the color term matrix B to zero. Residual systematics 
are clearly evident, in the form of deviations that ap- 
pear to increase with color. This is expected because we 
are comparing standard SDSS data to data taken with 
a different instrument than SDSS. We remedy this by 
measuring color terms for the instrument once using a 
standard stellar field observation, and then fixing these 
nonzero color terms while performing SLR. Color term 
issues are described in the following section. 

3.3. The Color Term Matrix 

The entries of the color term matrix B are zeroes or 
constants. For example, for the SDSS colors we use in 
i i4.2ll4.4[ c— {g — r,r — i,i — z), we adopt the color term 
matrix 

/bg -br \ 

B = 6, -b, ] . (3) 

V b,-bj 

This has a direct correspondence with traditional pho- 
tometric color term formulations, as shown in Appendix 

m 

Color terms, and the color term matrix, account for 
broad, differential instrument sensitivity differences that 
arise when data are acquired with different telescopes, 
CCDs, or filters than those used to generate the stan- 
dard catalog. For example, we see a clear systematic 
error between IMACS and SDSS in Figure [1] that varies 
monotonically with color. By estimating color terms with 
an observation of a standard star field and applying the 
color term matrix transformation during our SLR, we 
significantly improve the fit. 

We have also estimated color terms directly with SLR 
by letting the color terms be free in addition to k during 
the fit. This results in a divergent regression when using 
our weighted color residual GOF ( §3.6p . The reason is 
that the global minimum of this GOF occurs when all in- 
strumental data points collapse to a single point, which 
is allowed by divergent shears and rescalings from B, and 
divergent, compensatory color translations that put the 
singular instrumental data somewhere on the standard 
locus line. The best-fit color terms are extremely large, 
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violating our assumption of smallness. We have experi- 
mented with other GOFs with varying success, but none 
yields color terms with accuracy that rivals that of the 
traditional procedure. 

This is the only step in our real-time color calibra- 
tion procedure (i i3.8p that requires standard star obser- 
vations. If instrumental color terms are stable over, say, 
month- and year-long timescales, then color terms need 
only be estimated as infrequently. Without the require- 
ment of multiple standard field exposures per night, the 
observer maximizes the total exposure time on science 
fields. 

3.4. Color- Airmass and Higher Order Corrections 

Color-airmass terms take the same essential form as 
instrumental color terms. Appendix |A] gives an explicit 
example of this. Correspondingly, color-airmass terms 
can be estimated with the same color term procedure de- 
scribed in ^ 33.31 and §3.81 We have not yet implemented 
this as we assume these corrections to be small; this is an 
obvious future addition to our SLR formulation. In prin- 
ciple, corrections proportional to higher order powers of 
the color, airmass, and even the Galactic dust column 
can also be measured, but these will require larger, ded- 
icated programs in order to minimize error. As shown in 
i 34.2M.4( we achieve 1-2% level self-consistency with re- 
spect to zeropoints, airmass, and Galactic dust, without 
making these corrections. 

We now show how, using the same formalism and the 
same standard locus of §2.11 SLR can be made to out- 
put individual calibrations for each filter, producing cal- 
ibrated photometry instead of colors. 

3.5. Photometric Calibration with SLR Using 2MASS 

If the stellar locus is extended into other passbands 
that are already photometrically calibrated, then the in- 
strumental photometry can be directly calibrated using 
only the SLR methodology. 2MASS is an obvious choice 
for an external catalog, as it is full-sky and freely avail- 
able and the standard locus (i 32.ip bridges JHK with 
ugriz. 

The procedure requires the additional prior step of 
cross-correlating instrumental stellar catalogs against 
2MASS's calibrated data. This is easily done using the 
Gator web interface^. 

For illustration, say we have obtained the J-band pho- 
tometry of our instrumental griz stars. We construct a 
color that is a hybrid of instrumental optical magnitudes 
and calibrated 2MASS magnitudes, z — J, which extends 
the instrumental stellar locus toc = {g — r,r — i,i — z^z — 
J). SLR is executed in precisely the same way as for the 
optical data alone, using the standard optical-infrared 
hybrid locus line. The result is still /«, but the last en- 
try is Kzj = Uz + + Az because the 2MASS data are 
already calibrated: KzJ is the z-band photometric cali- 
bration. This can be done for any combination of instru- 
mental calibrated 2MASS data and instrumental optical 
or infrared data from a different passband. Appendix [Xl 
makes the mathematics of SLR photometric calibration 
explicit for this particular example. 

SLR photometric calibration is subject to the errors of 
2MASS, and is more precise when more stars are used 



in the fit. Given typical 2MASS errors of 5%, SLR- 
calibrated photometric zeropoints should typically be ac- 
curate to 5% or better. By the nature of SLR, the Galac- 
tic dust correction is inextricable from the photometric 
calibration when the instrumental stars are all behind 
the dust, so uncertainty from the dust must be added 
in quadrature. Finally, color terms carry with them ad- 
ditional uncertainty from the color term estimate and 
the colors that they multiply. If color terms are used, 
then the colors must be estimated independently from 
the the non-2MASS data, either simultaneously during 
the 2MASS fit or prior to it. 

2MASS matches will probably only occur for a sub- 
set of one's instrumental data. The reduced number of 
stars will generically degrade the Stellar Locus Regres- 
sion errors. Our optimized procedure is to first perform 
an instrumental-only SLR to estimate the colors using 
a maximal number of stars, and, if the photometry is 
needed at all, to perform a 2MASS SLR on the subset 
of matches to calibrate the photometry separately. This 
ensures that the color errors are minimal and not subject 
to the errors from 2MASS and from the reduced statis- 
tics. SLR never requires the subtraction of calibrated 
photometry to arrive at colors. 

3.6. Goodness of Fit 

During Stellar Locus Regression we optimize goodness- 
of-fit (GOF) statistics and estimate errors numerically. 
The GOF statistic we have adopted is the weighted, per- 
pendicular color-distance residual. 

The perpendicular color-distance is a hyper- 
dimensional distance in color space between an 
instrumental data point and the nearest point on the 
standard locus line. Our standard locus line is a collec- 
tion of closely-spaced data points, so we make numerical 
approximations to calculate the distance. As we make 
explicit below, we first compute all possible distances 
between the instrumental colors and the points on the 
standard locus line, then we find the minimum distance 
for each instrumental data point. Finally we sum the 
result. 

Consider instrumental color data points Cq,, where a 
indexes the instrumental color data vector of each star. 
We compute the vector distances between every Cq and 
every data point on our standard locus line, Cq/j, where 
Co/3 is also an array of vectors and /3 varies along the 
standard locus line. These distances are dafs = Cq, — Cq/j, 
where dap is a tensor of all possible distances between 
the instrumental colors and the points on the standard 
locus line. 

We weight the distances either by the number 1 or by 
the color measurement uncertainty in the direction of 
the line connecting each pair of points a-(3. Whether to 
weight by unity or by color uncertainties is at the dis- 
cretion of the user. We take the norms of each weighted 
distance vector, so our weighted distance measure is 

I da/3 I 



\d 



■a/3\ 



Here the dot product • is taken between the vector of 
data uncertainties <t„ and the unit vector that lies along 
the pairs of points da/j. The quantity cTq • d^p is the 
uncertainty projected along the line that connects each 

|http : //irsa . Ipac ■ caltech . edu/cgl-bln/Gator/nph-dd?catalog=f |Pfiae| of points. 
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Then, for each a, we select the one (3 = (3* such that 
|d^^. I — min(|<i*^|). This is the closest distance to the 
standard locus line for the instrumental data point a. 
We do this for all data points a, and sum the results as 

noting that there is one fi* for each a. This is an estimate 
of the weighted distance residual between all instrumen- 
tal data points and the standard locus line — a scalar. 

Our SLR implementation produces best-fit calibration 
parameters by varying k and adding it to the instrumen- 
tal colors, then recomputing T, and repeating this pro- 
cess until T is m inimized. We use the amoeba downhill 
simplex method (jNelder fc Meadlll965D . 

Our statistic was chosen for robustness. T would be 
equivalent to the statistic if the distances were squared 
and weighted by the inverse variance^, but the square 
puts undue weight on statistical outliers. Likewise, us- 
ing uniform weights of unity, the residual of square per- 
pendicular distances becomes equivalent to the so-called 
"total least-squares" statistic. This similarly gives undue 
weight to non-statistical outliers. We wish to apply SLR 
to any data at the telescope, on the fly, so robustness to 
outhers is critical. 

To insure doubly against outliers, we perform SLR 
twice. The first iteration gives rough estimates of best- 
fit calibration parameters, allowing us to excise stars 
we deem to lie too far away from standard locus line. 
We typically cut stars with color distance jdo,^. | > 

6|<Tq, • cZq,^. I aw ay from our stand ard locus line (a fig- 
ure informed bv lCovev et al.|[2007f ). and as well as those 
with Ida^g. I > 1 mag. We then perform SLR once more 
on the cleaned data. 

3.7. Uncertainty Estimation 

We estimate all uncertai nties on bes t-fit parameters 
with the bootstrap method (|Efronlll979f ). At each boot- 
strap iteration i we resample A^stars data points Ngtars 
times with replacement to obtain a bootstrap sample, 
and perform the GOF optimization described above on 
the new sample. This gives best fit parameters Ki, which 
we record. We repeat this A^boot 100 to 1000 times, 
obtaining A'boot-size arrays of each vector Ki. The mean 
of each bootstrap distribution is an estimate of the sam- 
ple mean of each parameter. The standard deviation of 
the bootstrap distribution is an estimate of the standard 
error on the mean. 

This approach means that to estimate uncertainties 
we repeat the entire stellar locus regression iVboot times. 
This can be time-prohibitive for large samples, but we 
argue it is not always necessary to repeat the error es- 
timation after it has been performed once. If different 
stellar images are taken from the same instrument under 
similar circumstances, and if the bootstrap distributions 
are roughly normal, then errors for different size data 
sets can be estimated with rescaling by yl^stars- This 
estimate may break down if those assumptions do not 
apply, but we argue V -/Vgtars rescaling can give errors to 

^ In this case our GOF is c l ose ly related to the 7-dimensional 
color distance of ICovev et al.l 112007) . The difference is that we 
weight by the errors in the direction of the line connecting the 
data and the nearest point on the standard locus line. 




Fig. 3. — SLR flow chart for calibrating colors. The hashed 
region denotes parts of the algorithm that are unique to SLR, while 
the non-shaded region shows steps that are more traditional. The 
dotted region denotes the color term estimation routine, which need 
only be performed once per detector. 

a factor of a few for Ngtars that are not wildly different. 
We now describe the full algorithm we have implemented 
to produce calibrated colors from flatfielded images. 

3.8. An Algorithm 

We have developed an optimized algorithm that pro- 
duces SLR-calibrated colors from flatfielded images, es- 
sentially in real time. It is schematically outlined in Fig- 
ure [3] The data undergo standard preprocessing, color 
terms are determined in the normal way from standard 
star frames once, and SLR fitting is done subsequently 
on all science frames. 

SLR requires single-epoch multiband observations of 
any number of science fields and one standard star 
field. All images are bias-subtracted, flat-fielded, WCS- 
registered, and Source Extracted. The resulting instru- 
mental catalogs from each passband are cross-correlated 
with one another to find the unique objects, and point- 
sources are identified, for example, using Source Extrac- 
tor's CLASS_STAR parameter. All stars are required to be 
unambiguously identifiable in all bands, and at all stages 
we only consider those with signal-to-noise > 10. 

Instrumental color terms are then determined once 
only from the standard star field, in the traditional way. 
For example, following the equations of Appendix [Al we 
estimate the g-band color term bg as the slope of the 
best-fit line to g — go vs. go — rg, and so forth for the rest 
of the passbands. The color terms are stored as the color 
term matrix B. 

With color terms in hand, real-time SLR calibration 
of any number of science fields can be undertaken. The 
color terms are first applied to the standard stellar lo- 
cus with the transform (1 -I- B), bringing the standard 
locus into standardized instrumental color space. In this 
space, the instrumental stellar colors from the flat-fielded 
science images are regressed to the standard locus, leav- 
ing only the color translation vector k free. This stage 
can be executed as soon as flat-fielded multiband images 
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are in hand. Flat-fielding and object cataloging can be 
performed at the telescope during observation, so Stel- 
lar Locus Regression is viable as a real-time calibration 
technique. 

The instrumental science catalog of both point-like and 
extended sources is calibrated by applying the inverse 
color transformation, 

Co = (l + B)-i(c-/^). (6) 

B was determined from the standard star observations, 
and K was obtained from the SLR. The transform (1+B) 
is invertible under normal circumstances of small color 
terms (see Appendix 1X1) . 

As outlined in SLR generates calibrated magni- 

tudes using 2MASS entries using exactly the color cali- 
bration process described above, with judicious choices 
of color vector entries. The additional step of match- 
ing one's instrumental objects to 2MASS lasts seconds 
to minutes for typical catalogs containing of order tens 
to hundreds of stars. 

We apply this algorithm to real data in a series of test, 
presented in the following section. 

4. FIRST TESTS 

We employed the algorithm of on existing, cal- 
ibrated SDSS data and on new griz data we have ac- 
quired on the 6.5 m Magellan telescopes using the IMACS 
(iDressler et"alll200l lOsip et al.l [20081 ) and LDSS3 (see 
lOsip et al.l 12001 instruments in imaging mode. We de- 
signed a series of tests to isolate the effect of Galactic 
extinction and airmass and to generally assess the repro- 
ducibility of colors using SLR. Finally, we used SLR in 
a measurement of redshifts of 11 galaxy clusters using 
colors alone. 
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Fig. 4. — Three samples of SDSS stars with different ranges 
of predicted Galactic extinction (low: Ag ~ 0.05 mag, intermedi- 
ate: Ag ~ 0.5 mag, and high: Ag ~ 1.2 mag). Top panels: color- 
color stellar loci, shown with the canonical reddening unit vector 
in three-dimensional color space. Bottom panels: histograms for 
predicted g-band Galactic extinction. 



TABLE 3 
Galactic Extinction Fields. 



4.1. Galactic Extinction 

The majority of SDSS stars at high Galactic latitude 
are expected to be behind the dust (^2]), so we designed 
a simple test of SLR to measure th e extinction directly. 
We a pplied SLR to ubercalibrated (Pad manabhan et al.l 
I2008D SDSS stars subject to varying degrees of predicted 
reddening, and compared the results to SFD expecta- 
tion. Ubercal SDSS data are cahbrated to the top of the 
atmosphere, in front of the Galactic dust. We therefore 
expect the best-fit color vector k to be equal to the red- 
dening vector, especially at high Galactic latitudes, and 
we expect the canonical Rv = 3.1 extinction law to hold. 

We began by querying photometric quality ubercal- 
ibrated SDSS point sources from the CasJobs web 
server^". Photometric quality ubercalibrated colors are 
accurate over a wide area of the sky, with 1.5% esti- 
mated uncertainty, so this is the ideal SDSS data set 
to perform our test on, in terms of size and quality. We 
queried stars in small ranges of SFD-predicted extinction 
and performed multiple such queries over a wide range 
of mean extinctions. The results were samples of hun- 
dreds of stars each in finite Galactic extinction bins over 
a wide range of SFD dust columns. We selected stars 
with signal-to-noise ratio > 10 in all bands to minimize 
uncertainty from poor color measurement. 

Figure 2] displays three of these data sets. Low 
{Ag ~ 0.05 mag), intermediate {Ag ^ 0.5 mag), and high 
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Mean Ag (mag) 


e (deg)=' 


b (deg)b 


0.048 


57.9 ± 1.0 


37.9 ±0.5 


0.538 


15.8 ±0.2 


29.0 ±0.2 


1.210 


11.4 ±0.1 


36.2 ±0.1 


0.090 


130.6 ±0.9 


-49.0 ±0.2 


0.268 


115.2 ± 0.4 


-48.7 ±0.1 


0.727 


149.8 ±0.4 


-46.6 ±0.1 


1.068 


56.6 ±33.9 


-37.5 ± 3.1 



^ Mean and standard deviation of Galactic lon- 
gitudes for stars in the sample.'' Mean and 
standard deviation of Galactic latitudes for 
stars in sample. 



{Ag ~ 1.2 mag) Galactic extinction ranges were queried, 
as shown in the figure. Histograms of the SFD-predicted 
extinction in the g-hand for each star show that these 
extinction distributions are localized and well separated. 
The gross effect of Galactic reddening is readily apparent 
in color-color diagrams (top panels). 

In total we acquired 7 such data sets, sampled from 
the northern and southern Galactic hemispheres, prob- 
ing mean Ag between 0.05 and 1.2 mag. By nature of 
the SDSS database organization and the available data 
itself, the data sets happen to come from fields that are 
localized to ~ Ideg^, except for one. The Galactic co- 
ordinates are given in Table [31 The last field in Table 
[21 shows a wide range of stellar coordinates because the 
CasJobs query, by chance, returned stars spread over a 
few disjointed regions of the sky. 

We performed Stellar Locus Regressions on each data 
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Fig. 5. — Best-fit color shifts k. vs. predicted Galactic extinc- 
tion in tlie g-band, and residuals. The solid lines are the SFD- 
predicted reddening corresponding to the g-band extinction given 
on the abscissa. Northern and southern Galactic hemisphere fields 
are indicated, and the itlOmmag region is shaded. 

set. We fixed the color term niatrix B to zero because 
this test involved comparing a standard stellar locus 
standardized to the SDSS photometric system, to col- 
ors generated by the SDSS instrument itself. After our 
restrictions on the data were enforced (see §3.6p , between 
380 to 490 stars from each data set were ultimately used 
in the fits. The best-fit color shifts k that resulted are di- 
rectly compared to the SFD-predicted reddening in Fig- 
ured The expectation, if our hypotheses hold, is that k 
will be equal to the reddening vector. 

The solid lines in Figure [5] show the 'predicted color- 
reddening values An — Am, n,ra £ {g,r,i,z\ correspond- 
ing to the predicted g-band extinction given on the ab- 
scissa, assuming the canonical Ry = 3.1 reddening law. 
The error bars for the abscissa are the standard deviation 
of predicted extinction values for the corresponding data 
set. Error bars for the ordinate are the standard error on 
the mean predicted reddening added in quadrature to the 
bootstrap errors of SLR and the 1.5% color uncertainty 
limit to which ubercalibrated colors are subject. 

The best-fit color shifts from SLR are consistent with 
Galactic reddening within the errors for these data, with 
residual RMS of 19mmag in the color g — r, ISmmag 
in r ~ i, and 13mmag in i ~ z. The maximum residual 
disagreement is 30 mmag in g — r and r — i, and 19 mmag 
i n i — z. 

llvezic et al.l ()2007l ) performed a similar analysis in 
SDSS Stripe 82 with their "stellar locus method," show- 
ing maximum sensitivity to the magnitude of SFD extinc- 
tion at roughly the 20 mmag level at high Galactic lat- 
itudes. Their measurements showed disagreement with 
SFD prediction at Galactic latitudes |6| ~ 40-50 deg and 
below (their §2.7.1). Our results probed latitudes in this 
range and lower, but in fields from different Stripes, so 
our results are not inconsistent with theirs. In other tests 
(i i4.2p . we observed Stripe 82 fields where their stellar lo- 
cus method failed, and we reproduced their results. Tak- 
ing this information together, we conclude that stellar 
locus methods can reproduce the Rv = 3.1 reddening 
law in both magnitude and direction, even through sig- 
nificant dust columns, so dust thickness alone is not a 
good indicator of SLR reliability. 

4.2. Color and Magnitude Reproducibility 

We tested the ability of SLR to reproduce colors and 
magnitudes using data acquired at different telescopes. 
We chose Stripe 82 fields because of the high quality. 



1% photometry available from llvezi(5 et al.l (|2007f ). We 
probed fields that contained among the highest densities 
of stars in the entire Stripe 82 in order to minimize statis- 
tical uncertainty. The stars in these regions, however, are 
known not to reproduce the SFD extin ction pre diction 
with stellar locus methods (§2.7.1 of ilvezic et al.i r2007). 
The possible reasons cited were an invalid extinction law 
Rv, a drift in mean metallicities, or dust that was both in 
front of and behind the stars. Nonetheless, we designed 
our experiment to control for the anomalous dust, as de- 
scribed below, so that we could study in isolation SLR's 
fundamental ability to reproduce colors and magnitudes 
when applied to data from different instruments. 

In 2008 we undertook observing programs at the 6.5 m 
Magellan telescopes using the IMACS instrument in 
imaging mode. We observed three, high stellar density 
Stripe 82 fields, which we label SI, S2, and S3. Our field 
SI is centered at {i,b) = (48.65151, -26.04729) deg, S2 
at (^,5) = (51.86849, -33.74441) deg, and S3 at {l,b) = 
(50.43107, -30.24483) deg. These fields are all at lower 
Galactic latitude \b\ than those of ! j4.11 except for the sec- 
ond field shown in Table[3 Each IMACS CCD is 4' x 8', 
arranged in two rows of four. Just as we have described 
in t ^3.81 we flat-fielded the data and produced instrumen- 
tal, multiband catalogs. We associated stars between 
passbands by requiring their positions to agree to better 
than 1" in radius. We th en found the c orresp onding pho- 
tometry in the catalog of llvezic et al.l ()2007f ) in the same 
way, selecting only those stars with more than 4 SDSS 
observations per passband and with signal-to-noise > 10. 

In order to compare the colors of stars directly, we 
treated the SDSS catalog data as pseudo-IMACS obser- 
vations: we extracted only the stars that also appeared 
in the 8 IMACS CCDs, in all griz bands. This was re- 
peated for the three fields. Because we wished to com- 
pare SLR magnitudes in addition to colors, we also found 
the 2MASS J-band photometry for the stars. Of order 
10-100 stars per CCD per field were matched across the 
IMACS, SDSS, and 2MASS data sets. 

Following the procedure outlined in §3.31 and §3.81 we 
measured the IMACS color terms from the SDSS stan- 
dard stars. We exploited the fact that we observed 3 
standard fields, and measured color terms independently 
in all of them, one CCD at a time, taking the average of 
the result. We applied the same mean color term correc- 
tion to all CCDs. This color term procedure gave smaller 
color and magnitude residuals than (1) using different 
mean color terms for each CCD and (2) using separate 
color terms for each field and CCD. 

First, we ran SLR to calibrate only the colors. We con- 
trolled for the anomalous stellar locus shifts by regressing 
both the calibrated SDSS data and our instrumental col- 
ors to the standard locus line. The best-fit n for the 
SDSS data showed color translations that disagreed with 
SFD p rediction by ^ 50- 100%, consistent with the find- 
ings of llvezic et~all ()2007t) . By applying SLR to both the 
instrumental and SDSS data sets, we effectively canceled 
the anomaly by subtraction. 

In a separate, subsequent step we calibrated the pho- 
tometry using SLR, fixing the optical-only color shifts to 
those measured in the first step, as described in §3. 51 and 
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^3.81 We solved for the SLR-standardized photometry as 
9q = 9 - Kgj - bg{gQ - ro) (7a) 
ro = r - Krj - br{ro ~ io) (7b) 
io ^ i - Ku - bi{iQ - zq) (7c) 
zq = z - K^j - b;,{io - zo) (7d) 
where the color terms 6„ were the measured values for 
the IMACS catalogs, and zero for the SDSS catalogs. 
For the colors on the right-hand side, we used the SLR 
calibrations from the first, optical-only color iteration. 
The 

i^gjj i^rj, stc, were those obtained from the second, 
optical-infrared magnitude iteration of SLR. 

Figures [6] and [7] show the resulting color and magni- 
tude residuals. The figures show a lack of pronounced 
systematic error in SLR colors as a function of effec- 
tive temperature, g — i. We assess that colors calibrated 
purely with SLR are reproducible between the SDSS and 
the IMACS instruments to (18,6, 5)mmag in the col- 
ors {g ~ r,r ~ i,i — z). Magnitudes obtained only with 
SLR, using 2MASS, are reproducible between IMACS 
and SDSS to (44, 25, 19, 9) mmag in the {g,r,i,z) pass- 
bands. These numbers do not include the Galactic dust 
uncertainty because we applied SLR to both data sets, 
but they do necessarily include the intrinsic airmass cor- 
rection that SLR makes. Our following test isolated com- 
pletely the effect of the atmosphere using IMACS data 
alone. 

4.3. The Atmosphere 

We again tested the reproducibility of SLR-calibrated 
colors, this time isolating and varying the effect of atmo- 
spheric extinction. An advantage to the SLR approach is 
that is does not assume any particular functional depen- 
dence of color variation with airmass (i 33.2p . We therefore 
expect precise airmass corrections from SLR. 

We isolated the effect of the atmosphere by perform- 
ing SLR on the same sets of stars observed through 
2 different airmasses using the same instrument, and 
compared directly the SLR-calibrated colors of the 
matched stars. We observed the IMACS fields de- 
scribed in §4.21 through multiple airmasses with the griz 
passbands. We included two additional fields, S5 at 
{e,b) = (122.46159, -63.21097) deg and CLl at {£,b) = 
(96.07818, -61.31761) deg— the latter of which we ob- 
served at three different airmasses. The field S5 is in 
Stripe 82, but did not contain enough matches between 
IMACS, SDSS, and 2MASS to include in the analysis of 
S21 The field CLl is not in the SDSS footprint. 

The data preprocessing was identical to that of ^4.21 
except here we applied the nonzero color term corrections 
to all IMACS data, and we did not cross-correlate with 
2MASS as we only wish to examine colors. 

Figure [S] shows the results. The left panel shows the 
weighted mean residuals of stellar colors per CCD per 
field, with bootstrap errors, plotted against the number 
of useful stars extracted from each CCD. The right hand 
panel shows the mean of each set of 8 data points per 
field, with error bars reflecting the RMS scatter, plot- 
ted vs. the airmass ratio of the two observations. CLl 
was observed at three airmasses, and the highest two 
airmasses are each paired with the lowest airmass in this 
figure. 

This shows that SLR colors are reproducible over a 
wide range of airmasses, with residual RMS scatter in 
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Fig. 6. — SLR color reproducibility, where SLR was applied 
both to SDSS and new IMACS observations. Color residuals are 
plotted vs. number of stars used per SLR, [left panels) and vs. true 
g — i color {right panels). Each data point on the left represents 
an SLR calibration done on one CCD, and residuals vs. g — i are 
averaged over the union of all stellar residuals across all fields and 
CCDs. The ±10 mmag region is shaded. 

{g — r,r — i,i — z) of (14, 8, 5) mmag when the number 
of stars used is > 30, and (20, 10, 8) mmag overall. Color 
errors do not correlate well with increased airmass. As 
is expected, there is a dependence of the color errors on 
the number of stars used per regression. The reduced 
of the fits are roughly in the range [1,3], with larger 
values for bluer bands. The larger iVstars data exhibit 
slightly larger , and smaller iVgtars data exhibit slightly 
smaller x^- This suggests that the bootstrap slightly 
overestimates errors at small A'stars? and slightly under- 
estimates errors at larger A^stars- This may be due to 
mildly non-Gaussian distributions. Overall, these results 
suggest that SLR systematic errors are not dominated 
by airmass issues for the fields we tested, but instead are 
dominated by the statistical noise inherent to the method 
as applied to the particular data we have used it on here. 

As our final test, we deployed SLR colors in a photo- 
metric redshift measurement, as distance estimation is 
one of the prime uses of accurately calibrated colors. 

4.4. Galaxy Cluster Redshifts 

We used SLR-calibrated colors to recover the red- 
shifts of galaxy clusters. As part of the 2008 observ- 
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Fig. 7. — SLR magnitude reproducibility, where SLR was ap- 
plied both to SDSS and new IMACS observations. A 2MASS SLR 
magnitude calibration was performed, again on both data sets as 
in Figure |6] with one data point per CCD. The ilOmmag region 
is shaded. 



ing program we mentioned in m.'2\ and i j4.31 we ob- 
serv ed known galaxy clus ters from the REFLEX cata- 



log (|Bohringer et al 



I2004D. primarily a subset that are 



also Abell clusters (jAbell et al.l |1989( ). These observa- 
tions were acquired with t he LDSS3 c amera at the Mag- 
ellan 6.5 m telescopes (see lOsip et al . 2008) in imaging 
mode. Because the redshift of the clusters are known, 
this serves as a test of our color calibrations within the 
context of a full scientific analysis. 

We estimated LDSS3 color terms using standard star 
fields. Individual SLR calibrations were then obtained 
directly from each of the single-epoch griz cluster fields, 
and the B and best-fit k applied to the galaxy col- 
ors without applying any further corrections for airmass 
or Galactic extinction. We then selected the brightest 
red-sequence cluster galaxies from color-magnitude di- 
agrams, spanning ^ 1-2 mag fainter than the brightest 
cluster galaxy (BCG) magnitude. This typically included 
of order 10 galaxies. We calculated the weighted mean 
of their g — i = {g — r) + {r — i) colors and estimated the 
cl uster redshift with the empirical color-to-redshift tables 
of lLopesI (|2007f) . For the low redshift range we probed, 
the g — i color of red cluster galaxies varies rapidly and 
monotonically with redshift and is an ideal redshift esti- 



FlG. 8. — Color residuals between SLR calibrations of fields 
imaged with IMACS at different airmasses. The ilOmmag region 
is shaded. Left hand plots: Residuals vs. number of stars used to 
perform each SLR. Right hand plots: Mean of residuals over the 8 
IMACS CCDs vs. airmass ratio between the two observations, and 
error bars are the RMS of residuals per field. 

mator. 

Table 0] summarizes the galaxy clusters we targeted, 
along with their Galactic coordinates and our results. 
The Abell systems 3693 and 3738 have been identified 
by previous authors as multi-redshifts systems. We were 
able to isolate the sub-systems using color-magnitude di- 
agrams and image inspection, and we estimated their 
redshifts independently. 

The results for the 11 clusters are plotted in Figure 
m The g — i cluster color errors SEgi were estimated 
as the standard error on the mean g — i color of the 
red sequence galaxies we selected, and cluster redshift 
errors were taken to be SEo; /3, taken fro m the slope of 
the color-redshift relation of lLopesI (|2007D . The reduced 
of residuals is 0.6, so this error slightly overestimates 
the scatter. The residual cluster redshift RMS is cr^ = 
0.007, or cr^/(l + z) = 0.6% for the range of redshifts 
we measured. We emphasize that this is the estimated 
redshift error not per red galaxy, but per cluster, each of 
which made use of an ensemble of red galaxies. 

The — 0.007 RMS error corresponds to g — i color 
errors of ~ 20 mmag, again estimating from the roughly 
linear relationship between redshift and g ~ i. This 
level of residual error is roughly in accord with the SLR 
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TABLE 4 

Summary of the Cluster Data and Results. 



Name References'' Redshift Airmass'^ Seeing" -Bg-i* N stars^ 6'^ 

Spec'' SLR,'^ (arcsec) (mag) (deg) 



ABELL 3668 


1,2 


0.1496 


0.156 ±0.011 


1.9 


2.0-3.0 


0.083 


58 


-32 


ABELL 3675 


1,2 


0.1383 


0.150 ±0.012 


1.9 


1.8-1.9 


0.079 


67 


-35 


RXC J2023.4-5535 


2 


0.232 


0.238 ±0.012 


1.7 


1.7-2.3 


0.102 


70 


-35 


ABELL 3693 


1,3,4 


0.091 


0.093 ±0.016 


1.5 


1.8-3.0 


0.060 


62 


-35 


ABELL 3693 


1,4,5 


0.123 


0.120 ±0.016 


1.5 


1.8-3.0 


0.060 


62 


-35 


ABELL 3739 


1,2 


0.1651 


0.177 ±0.010 


1.3 


1.6-2.1 


0.054 


63 


-42 


ABELL 3739 


1,6,7 


0.1786 


0.183 ±0.010 


1.3 


1.6-2.1 


0.054 


63 


-42 


ABELL 3740 


1,3 


0.1521 


0.141 ±0.007 


1.3 


1.4-2.7 


0.086 


39 


-42 


ABELL 3836 


1,3 


0.11 


0.120 ±0.007 


1.3 


1.4-2.7 


0.032 


32 


-51 


RXC J2218.6-3853 


8 


0.1379 


0.138 ±0.022 


1.5 


1.2-1.9 


0.026 


12 


-56 


ABELL 3866 


1,9 


0.1544 


0.156 ±0.008 


1.4 


1.1-1.9 


0.020 


22 


-57 


ABELL 3888 


1, 10 


0.152912 


0.159 ±0.016 


1.4 


1.3-1.9 


0.029 


28 


-59 


AM 2250-633 


9, 11 


0.2112 


0.204 ±0.014 


1.2 


1.1-1.4 


0.046 


26 


-49 



Refe renc es— ( 1) [A bcU ct al. (1989), (2) Bohrinscr ct al. (2004), (3) Strublc & Rood (1999), (4) Katscrt c t~aI1 
^MB), (5')IZaritsk v et al.. (.2006 '). (6') .Vogc s ct al. (1999,), (7),Schwopc ct al. (2000), (8) Crudda ce et al. (.2002) , (9) 
\dc Grandi ct al. (199^), ClOl IPimbblet et al. (2006), (ll HArp fc Madorc (1987).*" Spectroscopic redshift reported 
by references.'^ Photometric redshift based on SLR color corrections, with standard errors on the mean redshift 
estimates of the brightest cluster member galaxies. Airmass at which the images were obtained." Range of 
seeing in the multiband images.' SFD-predicted Galactic reddening in (g — i).^ Number of stars used to determine 
the SLR color correction.'' Galactic latitude of the cluster. 
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Fig. 9. — Photometric galaxy cluster redshifts from SLR- 
corrected colors vs. spectroscopic redshifts, and residuals. Right 
panel: Redshift residuals vs. expected g-band Galactic extinction. 
The red dashed line is the redshift error expected using colors that 
are not corrected for Galactic extinction. There is no evidence for 
a redshift residual correlated with Galactic extinction, indicating 
that SLR makes an appropriate correction. 

systematic errors due to Galactic dust and the atmo- 
sphere, estimated in ^AAl i )4.2[ and i j4.3[ and the residu- 
als are inconsistent with erroneous Galactic dereddening. 
If Galactic extinction were completely unaccounted for 
by our calibration method, redshifts would lie along the 
dashed line of the right-hand panel of Figure [HI assuming 
no other systematic effects. Because our redshift residual 
trend is not consistent with the dashed line, these results 
provide further evidence (see ^4.1[) that SLR accounts for 
Galactic extinction in these fields. 

5. DISCUSSION 



We have provided a new way to calibrate instrumen- 
tal photometry. It differs substantially from traditional 
techniques, and so a side- by-side comparison of SLR with 
traditional approaches is in order. This the first item of 
discussion. We then list in detail the unique advantages 
and limitations of SLR, and complete our discussion with 
tasks for the future. 

5.1. SLR vs. Traditional Photometric Calibration 
5.1.1. Traditional Photometric Calibration 

The traditional path to calibrated colors has been to 
first calibrate magnitudes, then to subtract them. Mag- 
nitudes are calibrated by tying instrumental flux to pho- 
tometric standard stars, which most often requires ad- 
ditional exposure time investment outside one's field of 
scientific interest. The resulting zeropoints must be in- 
terpolated in space and time to the science fields, but 
only after extinction by the atmosphere is modeled and 
measured — which necessitates further standard field ob- 
servations over a wide range of zenith angles. 

The analysis of these data to extract calibrated colors 
traditionally involves the steps of: 

1. Bias subtraction and flat-fielding; 

2. Extracting instrumental magnitudes for all images, 
in all bands; 

3. Using either observations of standard stars or pro- 
gram objects over a substantial span of zenith an- 
gles, in each band, to determine atmospheric ex- 
tinction coefficients, and correcting each observa- 
tion for atmospheric extinction; 

4. Assessing whether the atmospheric conditions 
during observation were photometric, and thus 
whether the zeropoint and atmospheric extinction 
interpolations are valid; 
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5. Correcting for the frame- by- frame difference be- 
tween PSF-fitting magnitudes and aperture pho- 
tometry; 

6. Determining and correcting for any differential sen- 
sitivity mismatch between the instrument and the 
desired photometric system, and applying these 
color term corrections to the photometry; 

7. Determining Galactic extinction in the direction of 
the observation from some external source such as 
SFD, and making an appropriate correction to all 
magnitudes; and 

8. Generating colors of sources of interest by subtract- 
ing these fully calibrated apparent magnitudes. 

The result of this procedure is a source catalog that 
contains magnitudes and derived colors. Besides being 
costly in telescope time, this procedure is also suffers 
from extra uncertainty from the explicit mathematical 
modeling of atmospheric attenuation. Moreover, the in- 
dividual magnitudes must each be calibrated to a better 
fractional precision than the desired color error. In the 
future, when high accuracy all-sky photometric catalogs 
of faint sources are readily available in all bands of in- 
terest, this process will be somewhat simplified. In the 
meantime we offer the SLR approach. SLR renders the 
preliminary magnitude calibration step unnecessary, in- 
stead allowing for the direct determination of colors of 
all objects of interest. SLR also provides a way to cal- 
ibrate the photometry by tying the magnitude scale to 
the all-sky 2MASS catalog, again without the need for 
extra standard star observations. 

5.1.2. SLR Photometric Calibration 

The SLR approach matches the distribution of stars in 
the instrumental color-color space to a standard stellar 
locus, allowing us to replace the traditional analysis path 
with a streamlined set of steps. 

1. For each previously uncharacterized instrument, we 
observe only one field that contains stars whose cal- 
ibrated magnitudes are known. This allows us to 
establish, in the traditional way, the instrumental 
color terms that arise from filter and detector dif- 
ferences. 

2. All subsequent images are bias subtracted and flat- 
fielded. 

3. Instrumental magnitudes are extracted for objects 
in each image. 

4. Objects are immediately cast in instrumental color- 
color space, and we determine the transforma- 
tions needed to bring the instrumental stellar locus 
into agreement with the universal, calibrated color 
properties of stars. 

5. This color transformation is then applied to all pho- 
tometry, producing calibrated colors for all mea- 
sured objects. 

6. If desired, calibrated photometric zeropoints are 
determined by bootstrapping the photometry to 



any stellar photometric catalog that overlaps the 
program fields. At the time of this writing, 2MASS 
is the obvious choice. 

The result is a catalog of colors, and optionally magni- 
tudes that are calibrated in a separate step and whose 
uncertainties do not factor into those of the colors. As 
a fundamentally different calibration technique, SLR has 
unique (in)sensitivities to astrophysical and instrumental 
effects. We break out the unique advantages and limita- 
tions these afford. 

5.2. Advantages 

SLR corrects for atmospheric extinction, even if time- 
variable. — Atmospheric extinction distorts the apparent 
colors of celestial sources, compared to what would be 
observed at the top of the atmosphere. To first approx- 
imation, this produces a translation of stars in color- 
color space. The atmospheric transmission function is 
expected to change as a function of time, zenith an- 
gle, and azimuthal angle as parcels of aer osols and wa- 
ter v apor travel and evolve in the sky (cf IStubbs et al.l 
I2007t ). Traditional airmass and zeropoint interpolations 
cannot account for this, unless specially designed sys- 
tems to monitor atmospheric conditions along the line of 
sight are deployed in parallel during observation. More- 
over, the ~ X^/^ airmass dependence of saturated water 
lines in the z-band is simply not modeled by most ob- 
servers (i j3.2p . SLR is unique because it naturally cor- 
rects for all of these additive effects at once with k, and 
is therefore insensitive to any additive mis-modeling of 
atmospheric extinction. We showed in !j4.2l and !j4.3l that 
SLR corrects for atmospheric attenuation in the z-band 
with sub-percent accuracy. Our SLR formalism also sup- 
ports color-airmass corrections, but we have shown that 
we already achieve high quality atmospheric corrections 
through a wide range of airmasses without these. 

SLR corrects for attenuation through clouds. — As long as 
the exposures are long enough to homogenize the ex- 
tinction due to clouds blowing across the images being 
analyzed, then the SLR technique will compensate for 
a common flux diminution across each frame. Integra- 
tion times that exceed 60 seconds typically satisfy this, 
even for images that span a degre e on the sky. The m ulti- 
epoch SDSS Stripe 82 analysis of llvezic et al.l (|2007| ) sup- 
ports the assertion that we can treat the effect of clouds 
as "grey" extinction, which SLR corrects for at once with 

K. 

SLR corrects for passband sensitivity differences, even across 
different instruments and telescopes. — An overall multi- 
plicative difference between the system throughput vs. 
wavelength of two different cameras or telescopes can be 
caused by filter transmission functions, detector quan- 
tum efficiency, and other instrumental properties. These 
are some reasons the nominal zeropoints between two fa- 
cilities are not equal. We have shown that SLR accounts 
for zeropoints in all tests we performed (SI)- Differen- 
tial variedion in system throughput vs. wavelength brings 
about disagreements in apparent magnitudes and colors 
that change with objects' color. The SLR technique takes 
these into account by applying traditional color terms, 
measured from infrequent standard star observations. In 
color-color space, these cause apparent colors to scale, 
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shear, and rotate. Our SLR methodology fully permits 
color term correction, and we have demonstrated this in 
2121 and implicitly in ^4.41 using two different cameras. 

SLR circumvents aperture corrections and photometric arti- 
facts from PSF variations. — Because the SLR technique 
maps observed, instrumental magnitudes onto a stan- 
dard stellar locus, even if photometry from a given im- 
age in a given passband has some systematic photom- 
etry error from PSF-related issues, the applied color 
transformation will correct for these artifacts. Differ- 
ences in aperture corrections between passba nds produce 
a common displacement in color-color space. Ilvezic et al.l 
|2007) showed that using the stellar locus to correct col- 
ors of galaxies works well for photometry obtained with 
the SDSS analysis pipeline. We have demonstrated in 
i j4.4l that applying SLR to Source Extractor instrumen- 
tal magnitudes produces photometric redshifts that are 
in excellent agreement with spectroscopic redshifts, with- 
out additional aperture corrections. 

SLR avoids the need to take more than a single calibration 
frame per instrument, per filter. — The universality of the 
stellar locus in effect allows the colors from a set of multi- 
band images to be self-calibrating. We find that taking 
only one set of multiband images of a calibration field 
is required, to determine the instrumental color terms. 
At all times thereafter, on timescales of months or years, 
and as long as there is a sufficient number of stars (> 10) 
in each image being analyzed, there is no need to obtain 
additional calibration frames. Our galaxy cluster obser- 
vations of ii4.4l used color terms measured only once in 
an SDSS standard star field. 

SLR corrects for Galactic extinction, independent of Av and 
Rv . — Whatever processes might bring about an overall 
observed displacement of the stellar locus, SLR will make 
the appropriate correction as long as it is locally well cor- 
related. This includes Galactic reddening effects when 
the stars are behind the dust. Because a single global 
color correction is applied to all the data, the cataloged 
objects should span a sufhciently small region on the sky 
to assure they have a common Galactic reddening. We 
showed in ^14. II that SLR can recover the canonical red- 
dening through very high and very low dust thicknesses 
to ~ 20mmag in color. We point out that our pho- 
tometric redshift analysis in i)4.4l made no explicit dust 
correction of any kind, and results were nonetheless con- 
sistent with an intrinsic Galactic dereddening by SLR. 
The anomalous extinction results of ^ 14. 21 are suggestive 
of possible deviations from the typically assumed redden- 
ing law — although drifts in mean metallicity cannot be 
ruled out as yet. 

SLR takes advantage of all stars in the images, producing a 
"democratic" color calibration. — For fields of view that are 
currently typical of astronomical imaging instruments, 
Galactic stars are approximately uniformly distributed 
across the image. SLR allows for the calibration of all 
objects in a field using the very stars appearing in that 
field. As long as there are no effects that vary across 
the field of view, the SLR color calibration technique 
inherently produces a homogenized calibration of colors 
across the frame. 



SLR uses a single common standard stellar locus, and thereby 
reduces sensitivity to systematic errors in photometric ze- 
ropoints. — The photometry from SDSS exhibits zero- 
point drifts at the few p ercent level across the sky 
(jPadmanabhan et al.l |200"8[) . If these cataloged magni- 
tudes are used to derive local colors, the colors inherit 
a position-dependent systematic drift that is the differ- 
ence between the underlying photometric zeropoint er- 
rors. SLR uses a single common standard stellar locus 
for all color corrections, so the derived colors should show 
reduced position-dependent systematic errors. The pre- 
cision of SLR-derived colors should be outstanding. The 
accuracy of SLR-derived colors depends upon the accu- 
racy of the standard stellar locus, the extent to which 
this locus is representative of the population being ob- 
served, and the nature of the Galactic dust extinction, 
as discussed in fj^l 

SLR is flexible. — The simple matching of the observed 
instrumental stellar locus to a standard one allows for a 
diversity of input catalogs. The input can be purely flat- 
fielded, instrumental photometry, and SLR will correct 
them for atmospheric and Galactic extinction, and ze- 
ropoints. SLR will equally allow input photometry that 
has already been corrected for some or all of these ef- 
fects, as previous authors have shown. We have applied 
SLR to pre-calibrated SDSS photometry (! j4.2p and to 
instrumental photometry from two different instruments 
(^ ^4.31 S3)- The method makes no distinction between 
the various sources of locus shifts in color-color space, so 
it can be applied to photometric catalogs generated up 
to any stage in the traditional analysis chain. One can 
even envision mixing catalogs by, for example, acquir- 
ing griz imagery of a field and performing SLR, then 
acquiring further 2;-band images of the same field and 
coupling the new instrumental z-band catalogs to the 
already-calibrated gri ones. The resulting k will produce 
the new z-band calibration, and should reproduce the gri 
calibrations. SLR can be optimized to give fast calibra- 
tions for a wide variety of optimized observing strate- 
gies. Of course, SLR will perform best when the prior 
reduction steps do not induce additional color scatter nor 
scalings, rotations, and shears. 

SLR is fast: mountaintop reductions can provide accurate 
colors, and hence photometric redshifts, on-the-fiy. — Opti- 
mizing the use of allocated telescope time is an important 
goal. By performing real-time SLR reductions of pho- 
tometric data, observers can determine when a desired 
signal to noise ratio is attained. A specific application of 
SLR by our group is the real-time determination of pho- 
tometric redshifts of large numbers of galaxy clusters — 
an application we illustrated in tj4.4l In upcoming runs 
we intend to use SLR to make adaptive adjustments of 
integration times while observing clusters of known po- 
sition but unknown redshift. 

5.3. Limitations 

As a different way of calibrating photometry, SLR also 
carries some unique limitations. The main limitation is 
that SLR performs better in obtaining highly accurate 
colors than magnitudes in each of the bands. By the 
same token, the colors that are produced by the SLR 
approach are only as good as the standard stellar locus 
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used as the calibration standard. Of particular concern 
is the need of identify a set of calibration stars that have 
suffered minimal Galactic extinction, and whose median 
metallicity reflects that of the stellar populations we ob- 
serve in practice. 

Another principal limitation is that SLR calibrations 
necessarily correct for Galactic extinction. To SLR, the 
dust correction is as natural as the atmospheric extinc- 
tion correction. This is because our locus is standardized 
to stars suffering minimal extinction, while stars we ob- 
serve in practice will almost always be behind the dust. 
It is, however, more common to calibrate photometry to 
the top of the atmosphere and only optionally apply an 
SFD Galactic extinction correction by hand, if needed. 
There is considerable uncertainty in our understanding of 
the dust, and this must be folded into our a priori color 
error estimates. In nearly all the fields we have studied, 
results suggest that SLR corrects extinction through a 
wide range of dust thicknesses with high accuracy. Users 
should take special care not to double-correct the extinc- 
tion when using SLR. 

Consequently, SLR will produce discrepant results if 
the sources of extinction vary significantly across the field 
of view. This introduces a trade-off between larger field 
size, which allows more stars to be included in the re- 
gression, and spatially constant extinction, which SLR 
corrects for best. We have applied SLR on small fields of 
view with as few as 7 useful stars, and the uncertainties 
in these cases are predominantly statistical. If SLR is to 
be applied to larger fields, one possible way to minimize 
extra scatter from Galactic extinction is to first apply the 
SFD dust correction for each stellar position, and then 
allow SLR to make the residual color calibration. For 
widely varying dust across a field, this should suppress 
the large scale gradients and thus unnecessary scatter in 
the stellar locus. 

A third source of concern, both for traditional methods 
and SLR, are systematic differences between point source 
and galaxy photometry. Ilvezic et"al] (|2007D and our own 
results demonstrate that applying a color correction de- 
rived from stars to galaxies does produce reliable photo- 
metric redshift estimates for red sequence galaxies. Some 
level of systematic error is nonetheless expected, because 
the spectral energy distributions (SEDs) of Galactic stars 
differ from the underlying SEDs of possibly redshifted 
galaxies of all types. Moreover, heavily extincted sources 
with different intrinsic SEDs will not suffer identical color 
shifts (McCalL2004 ). These are fundamental limitations 
of broadband photometry: as integrals over the prod- 
uct of a source's SED, atmospheric and Galactic dust 
transmission functions, and the instrumental sensitivity 
curve, broadband measurements entangle the true SEDs 
with intervening atten uation effects almost irreversibly 
(jStubbs fc Tonrvl[2006f l. This systematic error afflicts all 
photometry that uses stellar calibration standards, and 
is not unique to SLR. 

We finally note that, in the optical, SLR implicitly re- 
quires the 5-band. This is because the critical optical 
feature that makes our realization of SLR possible is the 
kink in the stellar locus at r—i ~ 0.7 — the stellar locus in 
the (r—i, i^z) plane is virtually featureless. The kink (1) 
gives the one-dimensional locus line a distinctive shape 
that uniquely locates it in color-color space, and (2) pro- 
vides a component of the locus that is nearly perpen- 



dicular to the redd ening from Galactic and atmospheric 
extinction (see also lStraizvs et al1ll998l ). A purely linear 
stellar locus line does not allow for a unique stellar locus 
regression solution, as an instrumental locus could slide 
freely along that line. Inspection of Figure [T] reveals the 
importance of having y-band data. 

5.4. Future Directions 

Standard locus refinements. — We consider the standard 
stellar locus used for this paper to be a starting point. 
SLR's performance rests on the standard locus being rep- 
resentative of the population that any given observation 
probes. It may be advantageous to develop different stan- 
dard loci for different expected populations — for exam- 
ple, one for low Galactic latitudes, or different standard 
loci for different exposure times. The various standard 
loci could be measured empirically, or the expected per- 
turbin g effects coul d be rn odeled analytically and applied 
to the lCovev et al.l ()2007l ) locus. 

A fuller tre a tment of Galactic extinction. — As pointed out 
bv lMcCaiH (|2004f ). extinction by Galactic dust produces 
a shift in color that depends on the underlying photon 
spectral energy distribution of the source. For the cur- 
rent implementation of SLR we have ignored this effect, 
which will produce shears, scalings, and rotations as well 
as a simple translation in color-color space. These are 
analogous to color-airmass effects. We also expect that, 
if large fields must be calibrated in one SLR pass, it may 
be advantageous to apply the SFD correction for each 
stellar position first, as mentioned in ^5.31 This may re- 
duce scatter in stellar loci constructed from large data 
sets with widely varying extinctions. 

SLR vs. Galactic extinction and metallicity. — This paper 
has shown that SLR accounts for the Galactic extinction 
in nearly all fields we have studied so far. The notable 
exceptions are th e Stripe 82 results in §4.2) which cor- 
roborate those of Ilvezic et al.l ()2007D . Interestingly, the 
Galactic latitudes of the Stripe 82 fields where stellar lo- 
cus methods fail are roughly equal to and higher than 
many of the fields we have used to recover the SFD ex- 
tinction (Snj and measure cluster redshifts (i j4.4p . This 
presents difficulty in placing a hard lower bound on |6| 
where SLR should be valid. While we expect SLR to per- 
form more reliably toward higher Galactic latitudes, we 
suspect SLR, by its nature, has the potential to make 
better Galactic extinction corrections than SFD. This 
warrants extended and in-depth studies comparing SLR 
color shifts jointly to SFD prediction and metallicity. 

Other passband combinations. — We have concentrated 
heavily on the grizJ bands, but the basic approach 
should be applicable to ot her photometric systems , such 
as u, HK, and UBVRI. iMacDonald et all (|2004D have 
shown that the stellar locus as probed by the Johnson 
bands has a kink feature that makes this possible. We 
will also be interested to see how the PanSTARRS y band 
at 1 /im will add to the accuracy of SLR. Another inter- 
esting path to pursue is the use of the SLR approach to 
transform data between photometric systems. 

Multi-frame SLR. — If we ignore color-airmass terms, 
then image stacks from individual frames taken at dif- 
ferent airmasses should also be amenable to SLR anal- 
ysis. A better approach, however, might be to extract 
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SLR colors from all possible independent permutations 
of the multiband images, and then average the result- 
ing colors. This is similar to the philosophy applied in 
the "N{N— 1)" approach to frame subtraction photome- 
try for supernova cosmology (|Barris et al.ll2005r ). wherein 
all possible pairs of images are subtracted. This seems 
especially worth pursuing for multi-epoch surveys like 
PanSTARRS and LSST, as it completely sidesteps the 
challenges of combining frames with different PSFs ob- 
tained at different airmasses. 

Color recalihration of SDSS and 2MASS?— If the stellar 
locus does in fact provide us with a uniform calibration 
source over most of the sky, and if SLR provides spe- 
cial insensitivities to anomalous atmospheric and Galac- 
tic extinction effects, then it is enticing to consider a 
joint recalibration of colors only from both the SDSS and 
2MASS surveys. A Bayesian approach using the unifor- 
mity of the s tellar locus, as was done in optical bands by 
llvezic et al.l (|2007t ). might improve the accuracy of both 
optical and IR colors in these catalogs. 

6. CONCLUSIONS 

We have developed and demonstrated a tech- 
nique that exploits the universality of the stel- 
lar locus to immediately obtain accurate colors 
from uncalibrated multiband data. For those 
who might wish to exploit the SLR approach, the 
core I DL tools we have developed are available at 
|http : // stellar-locus-regression. googlecode . com. 

Using archived photometry from the SDSS survey, we 
have demonstrated that SLR can produce results that 
agree with the commonly used SFD extinction map over 
a wide range of dust thicknesses. The performance of 
SLR in these cases appears to be limited by systematic 
zeropoint drifts in SDSS magnitudes. We also repro- 
duced the anomalous stellar locus reddening results of 
llvezic et al . (2007) in a few SDSS Stripe 82 fields. This 
puts into question (1) the validity of the Ry = 3.1 red- 
dening law in those fields, (2) whether the dust is both 
in front of and behind the stars, and (3) whether these 
fields have stars with spatially correlated deviations from 
expected median metallicity. 

Images of fields obtained through a wide range of air- 
masses were subjected to SLR analysis. We recovered 
the colors of the stars in these images with an uncertainty 
limited by the Poisson and flat-field errors in the photom- 
etry from each frame. Using SLR-only techniques and 
2MASS calibrated photometry, we additionally obtained 
z-band zeropoints in the fields good within 18 mmag. 

Finally, we have also presented photometric redshift 
results using SLR-derived colors only. We recovered the 
spectroscopic redshifts of 11 galaxy clusters at redshifts 
0.09 < z < 0.25 with (tJ(1 + z) = 0.6% RMS residual 
error. The redshift residuals also showed that the red 
sequence galaxy colors from SLR alone were consistent 
with an intrinsic Galactic extinction correction. 

The SLR technique, as implemented in our IDL code, 
can be used at the telescope in real-time to optimize the 



use of allocated telescope time. We regard SLR as a 
promising way to calibrate colors and magnitudes using 
fundamentally different physical assumptions, providing 
calibrations far faster than the traditional approach. 
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APPENDIX 

MOTIVATING THE COLOR TRANSFORMATION EQUATION 

Equation ([1]) represents photometric calibrations using a compact, filter-independent matrix notation because in the 
future we wish our procedure to be applicable to any set of filters. In this appendix we give a concrete example for 
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the griz SDSS filters in order to motivate and convey the meaning of the notation. 

Usual Photometric Calibration Equations 

Photometric cahbrations are normaUy modeled with simple additive and sometini es multiplicative terms. Following 
the southern SDSS standards literature (jSmith et al.ll2002l . [200l IClem et al.ll2007f) " as rough guides, we relate the 
instrumental magnitude {g, r, z, z) through the SDSS passbands to the true, extra-Galactic magnitude (goi ''Oi ^o) 



g = go + ag + Eg + Ag + bg{go - ro) + CgXg{go - ro) (Ala) 

r = Tq + Gr + E,. + Ar + br{ro — io) + CrXr{ro — io) (Alb) 

i = io + Oi + Ei + Ai + bt{io - zq) + c^Xi^io - zq) (Ale) 

Z ^ Zq + O;, + E^ + A^ + 62 (io - Zq) + C^X^{io - Zq) (Aid) 

J - Jo, (Ale) 



where in the last equality we take the 2MASS data to be calibrated already. Here, a„ is the zeropoint for the passband 
n; En is the atmospheric extinction, often modeled as fc,iX„ for some filter-dependent coefficient fc„ and airmass Xn 
through which exposure R is taken; and A„ is the estimated Galactic extinction. 6„ is the color coefficient, and c„ is 
the color-airmass cross term coefficient. The free parameters fc„, a„, 6„ and c„ can be estimated from the literature or 
measured using intermittent standard star observations, interpolating in airmass and in time to the science exposures. 
Galactic extinction is estimated using SFD. 



Corresponding Color Equations 
Colors are magnitude differences. Subtracting magnitudes between adjacent passbands using Equations 
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This is greatly simplified with matrix notation. We define the instrumental color vector as 




and likewise for the true color vector Cq; the additive color calibration vector as 



Eg — Er 



Ag Ar 



K. = 



the color term matrix as 
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and the color-airmass term matrix as 



fCgXg -CrXr 

CrXr —CiXi 

C,X, - CzXz 

V CzXz 



then we can write the color calibration equations as a single matrix equation, 

c = K+{l + B + T)co, 

For other relevant unpublished documents see |http : / / www-star . f nal . gov /Southern_ugriz/publications . htrnT] 



gives 

(A2a) 
(A2b) 
(A2c) 
(A2d) 
(A2e) 
(A2f) 
(A2g) 
(A2h) 



(A3) 



(A4) 



(A5) 



(A6) 



(A7) 
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where 1 is the identity matrix. Setting T = 0, this is equal to Equation ([1]). The inverse color transformation is 

co^{l + B + Ty\c~K). (A8) 

Whether the transformation matrices are invertible depends on the details of how first order color coefficients are 
chosen, but under normal circumstances it will be possible. Using the example we've just presented, (1 + B) is 
invertible because its determinant is (I + bg + br + bgbr){l + bi — bz) , which is not equal to zero for typical (small) 
values of color terms. 

This is how our matrix notation corresponds to typical photometric calibrations. The notation is trivially generalized 
to any filter set. The notation also motivates generalizations in choice of calibration parameters. For example, the 
color term matrix B can be made to be diagonal, (anti-)symmetric, or populated entirely by different nonzero entries, 
depending only the particular application. Furthermore, higher order terms can be considered by including higher 
powers of Cq to Equation (jA7p . 
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